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We report a measurement of the amplitude ratio rs of B° —y D°K*° and B° —y D°K*° decays 
with a model-independent Dalitz plot analysis using D —> Kgn + Tv~ decays. Using the full data 
sample of 772 x 10 6 BB pairs collected at the Y(4S) resonance with the Belle detector at KERB 
accelerator the upper limit is rs < 0.87 at the 68 % confidence level. This result is the first 
measurement of rs with a model-independent Dalitz analysis, and is consistent with results from 
other analyses. The value of rs indicates the sensitivity of the decay to 0 3 because the statistical 
uncertainty is proportional to 1 /rs- The rs result is obtained from observables (x±, y±) 


x- 

V- 

x+ 

y+ 


+0.4t 7 ;°±°° ± 0.0 

-0.6t?™ ±0.1 
+0.1±g;I±g;? dh 0.1 
-l-0.3±g;g±S:? ± 0.1, 


where x± = rs cos(5s ± 03 ), y± = rs sin(5g ± 03 ) and 03 (5s) are the weak (strong) phase difference 
between B° —y D°K*° and B° —y D°K*°. The first uncertainty is statistical, the second is the 
experimental systematic and the third is the systematic due to the uncertainties on a and Si 
parameters measured by CLEO. 


PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh, 14.40.Nd 
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INTRODUCTION 


Determination of the parameters of the standard model (SM) plays an important role in the search for new physics. 
In the SM, the Cabibbo-Kobayashi-Maskawa (CKM) matrix [lj gives a successful description of current measurements 
of CP violation. CP-violating parameters <j> i, </>2 and ^>3 are the three angles of a CKM unitarity triangle, of which 
</>3 = a,rg(—V u dVub*/V c dV c b*) is the least accurately determined. In the usual quark-phase convention, where large 
complex phases appear only in V u b and Vtd j2]> the measurement of (j )3 is equivalent to the extraction of the phase 
of V u b relative to the phases of the other CKM matrix elements except Vtd■ To date, <f >3 measurements have been 
made mainly with B meson decays into dV)rV) final states 0 - 12 |, which exploit the interference between the pK*)° 
and P>M° decaying into a common final state. In particular, Dalitz plot analyses of D decays in B± — > dV)rC)± 
processes are the most sensitive to <f >3 at e + e _ B-factories. The Dalitz plot analysis technique for the measurement 


of <563 was pro posed in Ref. [13, 14 ]. Belle reported the first </>3 measurement with model-independent Dalitz analysis 
technique Il5j . 

In this paper, we present the first measurement of the amplitude ratio of B° —> D°R*° and B° —>• D°R*° decays 
with a model-independent Dalitz plot analysis, which is essential for the determination of </>3 from this channel. We 
reconstruct B° —>• DR*° , with R*° —> R + n~ (charge conjugate processes are implied throughout the paper and 
R*° refers to KT*( 892 )°). Here the flavor of the B meson is identified by the kaon charge. Neutral D mesons are 
reconstructed in the R° s 7 r + 7r _ decay mode. The final states we reconstructed can be reached through b —>- c and 
b —> u processes with the diagrams shown in Fig. [l] 



Since the R* has a large natural width of 50 MeV, we consider interference between the signal B —>• DR* and 
B —>• DRtt where Rn arises from a non-resonant decay or through other kaonic resonances. In this analysis we use 
the variables r$, fc, and S 5 to parametrize the strong dynamics of the decay. These parameters are defined as [iff 


2 = T(B° —> D°R+n~) _ fdpAl_ u (p) 
rs T(B° — f D°R + n~) f dpA%_^ c (p) ’ 

ke iS s = I dpA b ^ c (p)A b ^ u (p)e iS ^ 

\f.I d P A b^c ( P ) I d P A b^u ip) 


where the integration is over the B° —> DR + ir~ Dalitz plot region corresponding to the R*° resonance. Here 
Ab^ c (Ab^ u ){p) are the magnitudes of the amplitudes for the b —> c(u) transitions, and 5(p) is the relative strong 
phase, where the variable p indicates the position within the DR + n~ Dalitz plot. If the B° decay can be considered 
as a DR*° two-body decay, r$ becomes the ratio of the amplitudes for b —> u and b —> c and k becomes 1 . According 
to a simulation study using a Dalitz model based on the measurements 


lit, 


the value of k is 0.95 ± 0.03 within 
the phase space of DR*° resonance. The value of rs is expected to be around 0.4, which naively corresponds to 
I V ub V* s | / | V c bV* s I, but also depends on strong interaction effects. 


THE MODEL-INDEPENDENT DALITZ PLOT ANALYSIS TECHNIQUE 

The amplitude of the B° —> DR*°, D Rgir + n~ decay is a superposition of the B° —> D°R*° and B° D°R*° 

amplitudes 


A_b(to+, m 2 _) = A + rse l< ' 5s+ ^ 3 ' > A, 


( 3 ) 
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where m 2 , and m 2 _ are the squared invariant masses of Kgir + and Kgir~ combinations, respectively, A = A(jn 2 + ,rn 2 _) 
is the amplitude of the D° —>• Kgir + n~ decay and A = A{m\,rriL) is the amplitude of the D° —> Kgir + n~ decay. In 
the case of CP conservation in the D decay ^4(m^_,TO?_) = A{rn 2 _,m\). The Dalitz plot density of the D decay from 
B° — DK*° is given by 

P B =| A b | 2 =| A + r s e* Sa+ +*)A | 2 = P + r|P + 2 kVpP(x+C + y+S ), (4) 

where P(rn\,m 2 _) =| A | 2 , P{m\,rrA.) =| A | 2 , and 

x + = r s cos(<5s + fa), y+ = r s sin(<5 s + fa). (5) 

Functions C = C(m+,rrvL) and S = S{mA,rn 2 _) are the cosine and sine of the strong-phase difference S B (rnTj_ . m 2 _) = 

arg A — arg A between the D° —> Kg n + n~ and D° —> Kg amplitudes. Here, we have used the definition of 
k given in Eq. ©■ The equations for the charge-conjugate mode B° —>• DK*° are obtained with the substitution 

fa •<-» —fa and A A; the corresponding parameters that depend on the B° decay amplitude are: 


x- = r s cos(<5 S - fa), y- = r s sin(<5 s - fa). 


( 6 ) 


Using B° and B° decays, one can obtain rs, fa and 5s separately. 

Up to this point, the description of the model-dependent and model-independent techniques is the same. The 
model-dependent analysis deals directly with the Dalitz plot density, and the functions C and S are obtained from a 
model based upon a fit to the D° 


—> KgTr + TT amplitude. In the model-independent approach EMI, the Dalitz plot 

The bin index i ranges from —TV" to TV" (excluding 
0); the exchange m 2 _ -o- mA corresponds to the exchange i -H- —i. The expected number of events in bin i of the 
Dalitz plot of the D meson from B° —y DK*° is 


is divided into 2TV" bins symmetric under the exchange mA -f-> m 2 , . 


TVf = h B K±i + r 2 gK Ti + 2k^K l K_ i (x±c i ± y±Si) 


(7) 


where stands for the number of B°(B°) meson decays, hs is a normalization constant and Ki is the number of 

events in the * th bin of the Kgir + Tr~ Dalitz plot of the D meson. The values of Ki are measured from a sample of flavor- 
tagged D° mesons obtained by reconstructing D* ± —» D-k± decays. The terms c,; and .s,; are the amplitude-weighted 
average of the functions C and S over the bin region: 


f Vi \A\\A\ CdV 
y/f Vi I 71 | 2 dDf Vi I A \ 2 dV 


( 8 ) 


Here D represents the Dalitz plot plane and T>i is the bin region over which the integration is performed. The terms 
Si are defined similarly with C substituted by S. The absence of CP violation in the D decay implies Cj = c-i and 
Si = —S—i. The values of Ci and s, can be measured using quantum correlated D pairs produced at charm-factory 
experiments operating at the threshold for DD pair production. The CLEO Collaboration has reported Ci and Sj 


values [ 20 U 21 I. Once the c,- and s,- measurements are included, the set of relations defined by Eq. © contain only 
three free parameters (x, y , and hs) for each B°, and can be solved using a maximum likelihood method to extract 
the values of fa, 6s and rs- We have neglected charm-mixing effects in D decays from both the B° —> DK*° process 
and in the quantum correlated DD production. 

In principle, the set of relations defined by Eq. © can be solved without external constraints on Cj and Si for TV" > 2. 
However, due to the small value of rs, there is very little sensitivity to the Ci and s t parameters in B° —>• DK*° decays, 
which results in a reduction in the precision on the other parameters. 


ANALYSIS PROCEDURE 

Equation 0 is the key relation used in the analysis, but it only holds if there is no background, no non-uniformity 
in the Dalitz plot acceptance and no crossfeed between bins. (Cross-feed is due to invariant-mass resolution and 
radiative corrections.) In this section we outline the procedures to account for the acceptance and crossfeed. 

First we discuss the effect of the variation of efficiency profile over the Dalitz plane. We note that Eq. © does 
not change under the transformation P —>• eP when the efficiency profile e(m^_,m?_) is symmetric: e(rrri,mA) = 
e(m 2 _, m+). The effect of non-uniform efficiency over the Dalitz plot cancels out when using a flavor-tagged D sample 
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with kinematic properties that are similar to the sample from the signal B decay. This approach allows for the removal 
of the systematic uncertainty associated with the possible inaccuracy of the detector acceptance description in the 
Monte Carlo (MC) simulation. With the efficiency taken into account (that is in general non-uniform across the bin 
region), the number of events detected is: 


N' = J p(V)e(V)dV. (9) 

Clearly, the efficiency does not factorize. One can use an efficiency averaged over the bin then correct for it in the 
analysis: 


Nj Jp(V)e(V)dV 
N t f p(V)dV 


( 10 ) 


The averaged efficiency Zi can be determined from MC. The assumption that the efficiency profile depends only on 
the D momentum is tested using MC simulation, and the remaining difference is treated as a systematic uncertainty. 
This correction cannot be performed in a completely model-independent way, since the correction terms include the 
amplitude variation inside the bin. Fortunately, calculations using the Belle D —> Kgi 


model [ 22 j show that 


this correction is negligible even for very large non-uniformity of the efficiency profile. 

There are two sources of cross-feed: momentum resolution and flavor misidentification. Momentum resolution leads 
to migration of events between the bins. In the binned approach, this effect can be corrected in a non-parametric 
way. The migration can be described by a linear transformation of the number of events in each bin 


N- = ^ a ik N k , 


( 11 ) 


where N k is the number of events that bin k would contain without the cross-feed, and N[ is the reconstructed number 
of events in bin i. The crossfeed matrix Oj k is nearly a unit matrix: on k -C 1 for i ^ k. The matrix is obtained from 
a signal MC simulation generated with the amplitude model reported in Ref. (22]. Most of the off-diagonal elements 
are null, only a few have values aii k < 0.04. In the case of a D —» /\]g 7 T + 7 r _ decay from a B , the cross-feed depends 
on the parameters x and y. However, this is a minor correction to an already small effect due to cross-feed; therefore 
it is neglected. 

The final effect to be considered are events in which the B flavor is misidentified. Double mis-identification in K*° 
reconstruction from K +/ rr~ , leads to migration of events between N+ •<->• iWj, due to assignment of the wrong flavor to 
the B candidate. If the fraction of double mis-identified events is /3, the number of events in each bin can be written 
as 


N* = N±+I3N*. ( 12 ) 

The value of /3 is obtained from MC data and is found to be (0.119 ± 0.007)%. Therefore, in this analysis the effect 
of flavor misidentification is neglected. 


Fit Procedure 

We can fit the data distribution in each bin separately, with the number of signal and background events as free 
parameters. The values of W found can be used in Eq. 0 to obtain the parameters (x±, y±). This is accomplished 
by minimizing a negative logarithmic likelihood of the form 

-2log C{x,y) = -2^1ogp((IVi) (x,y),Ni,a Ni ), (13) 

i 

where (Ay) is the expected number of events in the bin i obtained from Eq. 0. Here, Ay and a a are the observed 
number of events in data and the uncertainty on Ni, respectively. If the probability density function (PDF) p is 
Gaussian, this procedure is equivalent to a x 2 fit; however, the assumption of the Gaussian distribution may introduce 
a bias in the case of low statistics in certain bins. 

The procedure described above does not make any assumptions on the Dalitz distribution of the background events, 
since the fits in each bin are independent. Thus there is no associated systematic uncertainty. However, in the case 
of a small number of events and many background components this can be a limiting factor. Our approach is to 
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use the combined fit with a common likelihood for all bins. The background contribution has to be accounted for 
in the calculation of the values JV,. Statistically the most effective way of calculating the number of signal events 
is to perform, in each bin i of the Dalitz plot, an unbinned fit in the variables used to distinguish the signal from 
the background. In this analysis, we obtain the CP violating parameters (x ±, y±) from the data from a combined 
fit in bins. The relative numbers of background events in the bins of such a fit can be constrained externally from 
MC samples. In the combined fit, the expected numbers of events (Nj) as functions of (x, y) can be included in the 
likelihood. Thus the variables (x, y) become free parameters of the combined likelihood fit, and the assumption that 
the number of signal events has a Gaussian distribution is not needed. 


EVENT RECONSTRUCTION AND SELECTION 


This analysis is based on a data sample that contains 772 x 10 6 BB pairs, collected with the Belle detector 
at the KEKB asymmetric-energy e + e _ (3.5 on 8 GeV) collider [23| operating at the T(45) resonance. The Belle 
detector is a large-solid-angle magnetic spectrometer that consists of a silicon vertex detector, a 50-layer central 
drift chamber (CDC), an array of aerogel threshold Cherenkov counters (ACC), a barrel-like arrangement of time-of- 
flight scintillation counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux-return located outside of the coil is 
instrumented to detect K ° mesons and to identify muons. The detector is described in detail elsewhere [24]. 

We reconstructed B° -A DK*° events with K*° —> K + tt~ and D —> Kgn + TT~. The event selection, described 
below, is developed from studies of off-resonance data and MC simulated events. 

The Kg candidates are identified using the output of a neural network. Inputs to the network for a pair of 
oppositely-charged pions are the invariant mass, 20 kinematic parameters and particle identification (PID) information 
from the ACC, TOF and the ionization energy loss in the CDC. The Kg selection has a simulated purity of 92.2% 
and an efficiency of 75.1%. Charged kaon and pion candidates are identified using PID information. The efficiency 
is 85—95% and the probability of misidentification is 10—20% depending upon the momentum of the hadrons as 
obtained using dedicated data control samples. We reconstruct neutral D mesons by combining a Kg candidate with 
a pair of oppositely-charged pion candidates. We require that the invariant mass is within ±15 MeV/c 2 (±3er) of the 
nominal D° mass. K*° candidates are reconstructed from K + tt~ pairs. We require that the invariant mass is within 
±50 MeV/c 2 of the nominal K*° mass. We combine D and K*° candidates to form B° mesons. Candidate events 
are identified by the energy difference A E = — E b and the beam-constrained mass M^ c = ^ E\ — | JAp) | 2 , 

where £b is the beam energy and pi and Ei are the momenta and energies, respectively, of the B° meson decay 
center-of-mass (CM) frame. We select events with 5.2 GeV/c 2 < Mb c < 5.29 GeV/c 2 and 


products in the e + e 
-0.1 GeV < A E < 0.15 GeV. 

Among other B decays, the most serious background is from B° -A \K"*°t: + ] d + [K°n~] K ,- . This decay produces 
the same final state as the B° -A DK*° signal. To suppress this background, we exclude candidates for which the 
invariant mass of the K*°n + system is within ±4 MeV/c 2 of the nominal D + mass. This criteria leads to negligible 
contamination from B° -A [K* 0, k + \ d + [K 0 tt~] k *- and a relative loss in the signal efficiency of 0.6%. 

Large combinatorial background of true D° and random K + and tt~ combinations from the e + e _ -A cc process and 
other BB decays is reduced if D° candidates that are a decay product of D* + -a D°n + are eliminated. We use the 
mass difference AM between the [K { g’K + 'K~\g ) T: + and [iCg 7 r + 7 r _ ]£> systems. If AM > 0.15 GeV/c 2 for any additional 
7 r + candidate not used in the B candidate reconstruction, the event is retained. This requirement removes 19% of cc 
background and 11% of BB background according to MC simulation. The relative loss in signal efficiency is 5.5% 

In the rare case where there are multiple candidates in an event, the candidate with Mb c closest to its nominal 
value is chosen. The relative loss in signal efficiency is 0.8%. 

To discriminate the large combinatorial background dominated by the two-jet-like e + e _ — > qq continuum process, 
where q indicates u, d , s or c, a multivariate analysis is performed using the 12 variables introduced in the Table U) 
To effectively combine these 12 variables, we employ the NeuroBayes neural network package [27|. The NeuroBayes 
output is denoted as Cnb and lies within the range [— 1 , 1 ]. Events with Cnb ~ 1 are signal-like and events with 
Cnb ~ —1 are gg-like. The training of the neural network is performed using signal and qq MC samples. The Cnb 
distribution of signal events peaks at |Cnb| ~ 1 and is therefore difficult to represent with a simple analytic function. 
However, the transformed variable 


Cnb 


In 


Cnb - Cnb ,low 

Cnb ±igh Cnb 


(14) 






1 

Fisher discriminants based on modified Fox-Wolfram moments [251]. 

2 

The angle in the CM frame between the thrust axes of the B decay 
and that of remaining particles. 

3 

The signed difference of the vertices between the B candidate 
and the remaining charged tracks. 

4 

The distance of closest approach between the trajectories of the K* and D candidates. 

5 

The expected flavor dilution factor described in Ref. bdl. 

6 

The angle 9 between the B meson momentum direction and the beam axis in the CM frame. 

7 

The angle between the D and T(45) directions in the rest frame of the B candidate. 

8 

The projection on to the e + e _ beam direction of the sphericity vector 
with the largest eigenvalue. 

9 

The angle of the sphericity vector of the signal with respect to that of the remaining particles 

with the largest eigenvalue. 

10 

The angle of the sphericity vector of the signal with respect to the remaining particles 
with the second largest eigenvalue. 

11 

The angle of sphericity vector from the signal with respect to the remaining particles 

with the smallest eigenvalue. 

12 

The magnitude of the thrust of the particles not used to reconstruct the signal. 


TABLE I: 12 variables for qq suppression. 


where Cnb.Iow = —0.6 and C'NB.high = 0.9992, has a distribution that can be modelled by a Gaussian for signal as 
well as background. The events with Cnb < —0.6 are rejected. 

The number of signal events is obtained by fitting the three-dimensional distribution of variables M\, c , A E, and 
C^b using the extended maximum likelihood method. We prepare three-dimensional PDFs for each component by 
forming the product of one-dimensional PDFs for A E, Mb c and C^j B , since the correlations among the variables are 
found to be small. 

Backgrounds are divided into the following components: 

• Combinatorial background from qq background events. 

• BB background, in which the tracks forming the B° —> DK*° candidate come from decays of both B mesons in 
the event. The number of possible B decay combinations that contribute to this background is large, therefore 
both the Dalitz distribution and fit parameters distribution are quite smooth. In this analysis, BB backgrounds 
are further subdivided into two components: events reconstructed with a true D —>• Kg7r + 7r~ decay, referred to 
as D true, and those reconstructed with a combinatorial D candidate, referred to as D fake. 

• Peaking BB background, in which all tracks forming the B° —>• DK*° candidate come from the same B meson. 
This background has two types: events with one pion misidentihed as a kaon, such as Z?°[7r + 7r _ ] p o; and one pion 
misidentihed as a kaon and one pion is not reconstructed, such as Z? 0 [7r + 7T + 7r _ ] a +. 

The A E PDFs are parameterized by a double Gaussian for the signal, an exponential function for the D true 
BB background, an exponential function for the D fake BB background, a linear function for the qq background, a 
double Gaussian for the D°p° background, and a Gaussian for the D°a background. The M^ c PDFs are a Gaussian 
for signal, a Crystal Ball function for D true BB background, an Argus function for D fake BB background, an 
Argus function for qq background, a sum of Gaussian and Argus function for D°p° background and a Gaussian for 
D°ai background. The C^f B PDF are a sum of Gaussian and bifurcated Gaussian for each component. The shape 
parameters of the PDFs are fixed from MC samples. 

The Dalitz plot distributions of the background components are discussed in the next section. The numbers of 
events in each bin can be free parameters in the fit, thus there will be no uncertainty due to the modeling of the 
background distribution over the Dalitz plot in such an approach. This procedure is justified for background that is 
either well separated from the signal (such as peaking BB background), or is constrained by a much larger number 
of events than the signal (such as the qq background). The results of the fit over the full Dalitz plot are shown in 
Figs. [2j We obtain a total of 44.2l^ i signal events. The statistical significance is 2.8er. 
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FIG. 2: Projection of the fit to real data using the full Dalitz plot. Left: A E distribution with Mb c > 5.270 GeV/c 2 and 
Cnb > 2 requirements. Middle: C^b distribution with | A E |< 0.030 GeV and Mbc > 5.270 GeV/c 2 requirements. Right: 
Mbc distribution with | A E |< 0.030 GeV and C^b > 2 requirements. Histograms show the fitted signal and background 
contributions, (red is signal, yellow is D°af, green is D°p°, blue is D fake BB, light blue is D true BB and magenta is qq) 
and points with error bars are the data. 


Parameter 


y- 

corr.(x, y-) 

+0.29 ±0.32 

-0.33 ±0.41 
+7.0% 

x+ 

y+ 

corr.(x+, y+) 

+0.07 + 0.42 
+0.05 ± 0.45 
-7.5% 


TABLE II: ( x, y) parameters and their statistical correlations from combined fit of the B° —>• DK*° sample. The error is 
statistical uncertainty. The values and errors are obtained from likelihood distribution. 


DATA FIT IN BINS 

A combined fit is performed to obtain the B° -+ DK*° yield in each bin. The combined fit constrains the amount 
of the D true BB background in bins from the ratio of D° (Ki ) and D° {K_i) from the generic MC, the amount 
of the D fake BB, qq, D°p° and D°a+ backgrounds in bins from the MC, and takes the (x±, y±) variables as free 
parameters. Fits to B° and B° data are performed separately. The plots illustrating the combined fit results are 
given in the Appendix. In this lit, the additional free parameters are the total yields of D true BB, D fake BB, 
qq and peaking BB backgrounds for all over the Dalitz plane. The values of (x, y) parameters and their statistical 
correlations obtained from the combined fit for signal sample are given in Table [TT] 


SYSTEMATIC UNCERTAINTY 

The systematic uncertainties on (x, y) are obtained by taking variations from the default procedure under differing 
assumptions. Most systematic uncertainties are negligible compared to the statistical uncertainty. Therefore, numer¬ 
ical values are not given unless they are greater than or equal to 0.1 for any contribution. There is an uncertainty due 
to the Dalitz plot efficiency variation because of the difference in average efficiency over each bin for the flavor-tagged 
and B° -+ DK*° samples. A maximum difference of 1.5 % is obtained in a MC study. The uncertainty is taken as 
the maximum of two quantities: 

• the root mean square of x and y from smearing the numbers of events in the flavor-tagged sample Ki by 1.5%, 
or 

• the bias on x and y between the fits with and without efficiency correction for Ki obtained from signal MC. 
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The uncertainty due to crossfeed of events between bins is estimated by taking the bias between the fits with and 
without the correction. The uncertainties due to the fixed parametrization of the signal and background PDFs are 
estimated by varying them by ±ler. The uncertainty due to the C^j B PDF distributions for BB are estimated by 
replacing them with the signal C^ B PDF. The uncertainty due to D true and D fake BB fraction of yield is estimated 
by varying then from 0 to 1. The uncertainty due to PDFs shape is (Ax_, Ay_, Ax + , A y + ) = (lo;?,l°'J l >-oi )• 
The uncertainty arising from the finite sample of flavor-tagged D —> Kgirn decays is evaluated by varying the value 
of Ki within their statistical uncertainties. The uncertainty due to the uncertainty on k in Eq. m is evaluated by 
varying the value of k within its error & The uncertainty due to the limited precision of Ci and Si parameters is 
obtained by smearing the Ci and Si values within their total errors and repeating the fits for the same experimental 
data. The uncertainty due to Ci , Si is (Ax_, Ay_, Ax + , Ay + ) = (±0.0, ±0.1, ±0.1, ±0.1). 

The total systematic uncertainty is (Ax_, Ay_, Ax + , Ay + ) = (!?'?, ±0.1, ±0.1, ±0.1). The systematic uncertainty 

without r s is ('+ 0 0 + 0 - 1 +0.0 +0.0 \ 
witnoui c,, s t is t_Q i,_o.o >—o.i i—o.i )• 


RESULT FOR (x, y) AND rs 

We use the frequentist approach with the Feldman-Cousins ordering [2!§] to obtain the physical parameters 
p = (0.3, rg. <5s) (or true parameters p = Ztrue = (x-, y_, x+, y+)) from the measured parameters 3; = z mea s = 
(x_, z/_, x + , y + ) taken from the likelihood distribution. In essence, the confidence level a for a set of physical param¬ 
eters /i is calculated as 


, , /© (M) P( Z \^) dz 

a ^ focP( z I P"> dz 


(15) 


where p(z \ p) is the probability density to obtain the measurement result z given the set of true parameters p. The 
integration domain T>(fx) is given by the likelihood ratio (Feldman-Cousins) ordering: 


P{z I k) 

p{z I Pbest{z)) 


p(zq I P) 

p(z 0 | /Tbest(-Zo))’ 


(16) 


where /ibest (z) is p, that maximizes p(z \ p) for the given z, and zo is the result of the data fit. This PDF is taken 
from MC pseudo-experiments. 

Systematic errors in p are obtained by varying the measured parameters z within their systematic errors assuming 
a Gaussian distribution. In this calculation we assume that the systematic errors are uncorrelated between the B° 
and B° samples. 

As a result of this procedure, we obtain the confidence levels (C.L.) for ( x , y) and the physical parameter rg. The 
C.L. contours on (x, y) are shown in Fig. [3] The likelihood profile for rg is shown in Fig. [3] The final results are: 


X- 

= ±0.41 ± 0.0 

(17) 

y- 

= —0.61?™, ±0.1 

(18) 

x+ 

= ±0.l±8;I±g;S ± 0.1 

(19) 

y+ 

= ±0.31°™ ± 0 . 1 , 

(20) 

rg 

< 0.87 at 68% C.L.. 

(21) 


CONCLUSION 

We report the result of a measurement of the amplitude ratio rg using a model-independent Dalitz plot analysis 
of D — y KgTT + Tr~ decay in the process B° —> DK*°. The value of rg indicates the sensitivity of the decay to (f >3 
because the statistical uncertainty is proportional to 1/rg. The measurement was performed with the full data sample 
of 711 fb -1 (772 x 10 6 BB pairs) collected by the Belle detector at the T(4S) resonance. Model independence is 
achieved by binning the Dalitz plot of the D —>• K^-k + t:~ decay and using the strong-phase coefficients for bins 
measured by the CLEO experiment |2l|. We obtain the value rg < 0.87 at 68 % C.L. 
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FIG. 3: C.L. on ( xy~) (blue) and (x+, y +) (red). Dots show the most probable ( x, y ) values, lines show 68 % contours. The 
fluctuations come from statistics of pseudo-experiments and C.L. step used. 




This analysis is the first using the model-independent Dalitz technique on neutral B mesons. This measurement 
has resulted in lower statistical precision than the model-dependent measurement from BaBar with the B° —>• DK*° 
mode Q despite the larger data sample due to the smaller B° — > DK*° signal observed. The result is consistent 
with the most precise rg measurement reported by the LHCb Collaboration Q of rg = 0.240^^048 which uses 
B° —> [K + K~,K ± Tr^,n + n~]r)K* 0 decays. 
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FIG. 5: Projections of the combined fit of the B° —» DK*° sample on A E for each Dailtz plot bin, with the Mhc > 5270 MeV/c 2 
and Cnb > 2. requirements. The fill styles for the signal and background components are the same as in Fig. [2] 


APPENDIX 

The results of the combined fit to B° —»• DK*° and B° —»• DK*° samples separately for each bin of the Dalitz plot 
are shown in Figs. [5] and [6l respectively. The plots show the projections of the data and the fitting model on the AE 
variable, with the additional requirements Mhc > 5270 MeV/c 2 and C^g > 2. 
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